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Abstract 



^ijT I We derive a new parton-like formula, which establishes a simple connection 

Q^ ■ between the electroweak decay rate r{B ^ ^sj) and the rate of a free 6-quark 

(— I , decay. The main features of our approach are the treatment of the 6-quark as 

CL|| an on-mass-shell particle and the inclusion of the effects arising from the 6-quark 

O ■ transverse motion in the i?-meson. Using various 5-quark light-front (LF) distri- 

O . bution functions, both phenomenological one and the ones derived from current 

constituent quark models, and neglecting perturbative corrections we compute the 

photon energy spectra and the moments of the shape function. It is shown that 

S^ . the parton-like approach is fully consistent with the Heavy Quark Effective Theory 

(HQET) provided the 6-quark constituent mass is redefined in the way similar to 

that used in HQET to define the pole mass of the b quark. In this way the correction 

to first order in l/rrih can be eliminated from the total width in agreement with the 

general statement of HQET. We have also found that the photon energy spectra 

calculated in the LF approach agree well with the ones obtained in the ACM model, 

provided the same distribution function is used as input in both cases. In spite of 

the simplicity of the model our results show a fair good agreement both with the 

HQET predictions and available experimental data. 
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Studying of the photon spectrum in the weak radiative B — > Xgj is important for under- 
standing how precisely the total rate can be predicted in the presence of an experimental 
cut on the photon energy E. Experimentally, because of backgrounds, only the high en- 
ergy part of the photon spectrum can be detected with the present experimental cut 
E > 2.1 GeV at CLEO |I|]. Gross features of the spectrum such as the average photon 
energy may be used to measure the fundamental HQET parameters [|], 0. Since the 
6-quark is heavy compared to the QCD scale, the inclusive B -^ Xsj decay rate can 
be calculated in a systematic QCD-based expansion [|]. However, near the end point 
an important non-perurbative effect due to the soft interactions of the 6-quark with the 
light constituents has to be included. This so-called "Fermi motion" can be included in 
the heavy quark expansion by resumming an infinite set of leading-twist contributions 
into a shape function F{x) 0, where a scaling variable is defined as x = {2E — mb)/A 
with A = Mb — rrib- This is quite analogous to what happens for the structure function 
in deep-inelastic scattering in the region where the Bjorken variable xs — > 1. A model 
independent determination of the shape function is not available at the present time, 
however it may be possible to address this issue using lattice QCD 0. Ansatz for the 
shape function constrained by the information on its few first moments has been recently 
used in Ref. including the full NLO perturbative QCD corrections. 

As to phenomenological analyses of the photon spectra up to now they have been solely 
based on the ACM model 0,0. This model treats the heavy hadron as a bound state 
of the heavy quark and a spectator, with a certain momentum distribution. A light-front 
(LP) approach to consideration of the inclusive semileptonic transitions was suggested in 
Refs. |l^-[l^ and has been recently refined in |jl3[. The corresponding ansatz of Ref. [|13 



reduces to a specific choice of the primordial LP distribution function |V'(^,pi)P, which 
represents the probability to find the b quark carrying a LP fraction ^ and a transfer 
momentum squared Px = IP-lP- As a result, a new parton-like formula for the inclusive 
semileptonic 6 — »■ c, m width has been derived |jl3|], which is similar to the one obtained 



by Bjorken et al. |TJ] in case of infinitely heavy h and c quarks. 

Some of the dinamical features of this model get obscured by the integration over the 
lepton energy. They are seen in a cleaner way in the spectrum of the photons in the 



radiative B -^ X^^ transitions. In this paper, we extend the work of Ref. |]T3 to compute 



the non-perturbative corrections to the photon spectrum and the B — > Xs7 inclusive 
rate0. We strive to implement the binding and the 5-meson wave function effects on the 
photon energy and the invariant mass distributions of the hadrons recoiling against the 
photon. We will also study the comparison between the photon spectra dV/dE calculated 
in the LP and ACM approaches and will show that the discrepancy between the two is 
very small numerically. 

To begin with, we briefly discuss a derivation of the inclusive photon energy spectrum 
for the decays B — > Xsj in the context of the LP approach of Ref. ||13[. Similar to the 



ACM model the LP quark model treats the beauty meson B as consisting of the heavy b 
quark plus a spectator quark. Both quarks have fixed masses, rrif, and rrisp, though. This 
is at variance with the ACM model, that has been introduced in order to avoid the notion 
of the heavy quark mass at all. 

The 6 — >■ S7 transition at the fundamental level is generated by electroweak penguins 
I6| , |T^ . In the leading logarithmic approximation the decay b ^ s'y is described by the 



*A preview of this work can be found in 



effective Lagrangian 

AGf 
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where 



Ceff = -^=-VtbVtlc7{iJ,)mb{n)07, 



Or = ^^u,a^,{l + ^,)ubF>'''. (2) 

In Eqs. (|l]),(@) Gp is the Fermi constant, Vij are the elements of the CKM matrix, fhh{^Ji) 
and c-j{ix) are running 6-quark mass and the Wilson coefficient, respectively, evaluated 
at a subtraction point yU, and F^^ is the electromagnetic field strength tensor. The 
strange quark mass will be neglected throughout this paper; it only enters the final results 
quadratically as ml/ ml. For simplicity we neglect contributions of other operators in the 
effective Lagrangian which appear at next-to-leading order. 



To calculate the decay rate B — > Xg^ we use the approach of Ref. [|T^ which is based 
on the hypothesis of quark-hadron duality. This hypothesis assumes that the sum over all 
possible strange final states Xg can be modeled by the decay width of an on-shell h quark 
into on-shell c-quark weighted with the 6-quark distribution function /(^,pi). Going 
through the intermediate steps (for the detail see Ref. |T^) we obtain for the partial 
decay rate of the inclusive decay B — >• Xg-^ 

dT = 4xoToJ^-^f{^,pl)dT, (3) 

where 

To = ^,\Vt,V,s\'cUm,)m,{f,Yml (4) 

is the contribution of the matrix element of O7 to the b —>■ s'j decay rate, and xo = irib/MB- 

In Eq. @) the factor m^ comes from the two-body phase space. In HQET irib is usually 

associated with the h quark pole mass mib^poie- In our phenomeno logical consideration we 

associate mf, with the constituent mass of the 6-quark, see below. The normalization of 

fi^,pl) reads 

1 00 

7rJd^Jdklf{tkl) = l. (5) 





The factor 1/^ in Eq. (^ comes from the normalization of the B —^ bd vertex |T^. The 
phase space factor dr is given by 

dr = S[{pb - qf]EdE. (6) 

We choose the z-axis parallel to the 3-vector q, so that g+ = 2E, g_ = 0, where q± = 
qo ± Qz, then dr takes the form 



2 pi + mi. 



dr = 5(mi - ^^ '''"' )EdE. (7) 



In the first approximation we neglect p\ in the argument of the 5-function. Then, intro- 
ducing the scaling variable y = 2E/Mb, the photon spectrum dT/dE in B —>■ Xg'j takes 
the simple form 

1 dT{B -. X,7) 



Fo dy 



Ririy), 



where 

RLF{y) = -yfiy), (9) 

Xq 

_ oo 

and /(^) = T^ j dp\f{^,p\). Equivalently, one can write the spectrum in the standard 
QCD form |T§ 

Therefore the specific choice adopted for f{S,) corresponds to a particular form of the 
QCD shape function F{x) = {A/MBXQ)yf{y), where x = {y — Xo)/(l — a;o)-Q 

Following Ref. [|13[ we now account for the transverse motion of the 6-quark in Eq. 



(0) to find a bit more complicated expression for Rip{y) ||U 

RLF{y) = 2mi,7T f' f{^,pf)d^, (11) 

where the integration limits follow from the condition p*^ > 0, with p*^ = p*_^{^,E) = 
ml{^/y — 1). In this case the shape of the spectrum is obtained by direct integration of 
the distribution function. We shall see below that the difference between Ripiy) given 
by Eqs. (^ and (0) is very small numerically. 

In the free quark approximation, f{^,p'±) = 6{^ — C,o)S{p']_), the total inclusive width 
r{B —>■ Xgj) is the same as the radiative 6 — > S7 width Fq, and the spectrum of photons 
is a monochromatic line: 

-d^ = ^oSiE--). (12) 

The delta-function of Eq. (|l^) is transformed into a peak of a finite width due to the heavy 
quark motion. This effect is solely responsible for the filling in the windows between m;,/2 
and the kinematical boundary in the B meson decay, E^ax = Mb/2. 0. The expressions 
in Eqs. @, (|ll]) exhibit a pronounced peak which is rather assymetric. It is gratifying 
feature of the LF model, since it is in qualitative accord both with findings in QCD and 
experimental data. The perturbative corrections arising from gluon Bremsstrahllung and 
one-loop effects [^] also lead to a nontrivial photon spectrum at the partonic level. Since 
our primary object here is to discuss non-perturbative effects due to the Fermi motion, 
we will implicitly ignore perturbative gluon emission throughout our analysis. In this case 
the parton matrix element squared is a constant and can be taken out of the integral in 
Eq. d). 

Since we do not have an explicit representation for the B-meson Fock expansion in 
QCD, we shall proceed by making an ansatz for the momentum space structure of the 
wave function. This is model dependent enterprise but has its close equivalent in studies of 
B -^ Xg'j using the ACM model. In what follows, we will adopt both a phenomenological 
LF wave function and the LF functions corresponding to the various equal time (ET) 
quark model wave functions. As to the phenomenological ansatz, we use a model first 



^If one uses the Infinite Momentum Frame prescription ||l^, [|llj pt, — £,Pb, i-e. "T-b(0 = C^-^^sj then 
from Eqs. (|),(§) one easily derives dV ex E^f{y)dE [|o|. 

^The true endpoint is actually located at [Mg — {rriK + mTr)'^]/2MB ~ 2.60 GeV, i.e. slightly below 
Mb/2 w 2.64 GeV. 



written in Ref. |2^, and also employed in Refs. [|TT|,|T^ to implement the bound state 
effects in S-meson decays. It is written in the Lorentz-invariant form 



V'(IpI) =^exp{--VBVsp) 



L TTlsp 



(13) 



where vb and Vgp are the 4-velocities of the 6-quark and the quark-spectator, respectively, 
and Ep = a/IpP + 'mlp is the energy of the spectator. We shall use the normalization con- 

dition J p'^dp4)'^{\p\)/2£p = 1, in which case A/"^ = 2A/(m^ i^i(A)), where -ft'i(A) is the 



McDonald function. The function $(p^) = ip'^{\p\)/2ep represents a momentum distri- 
bution of the spectator quark in the B meson rest frame. We convert from ET to LF 
momenta by leaving the transverse momenta unchanged and letting 



1 



Pi 



(pf - Pi ) 



1. + _ Pl± + m'i 
r^' Pt ^ 



(14) 



for both the 6-quark {i = b) and the quark-spectator {i = sp). The longitudinal LF 
momentum fractions ^j are defined as ^sp = pfp/Ps, ^b = pt / ^bi "with p'l +p^ 



P+. In 



the 5-meson rest frame P\ 



B 



normalized according to Eq. 



Mb- Then for the distribution function f{^,p\) (^ = C,b) 
one obtains 



/(e,Pi) = ^ 



in 



1-e 



^f^ 



^7r-(l-0 



exp 




(15) 



where ^o = t^sp/Mb- The function ([l5|) is sharply peaked ai p]_ = Q, ^ = ^o- In what 
follows, we shall refer to the LF wave function of Eq. (|1^) as the case A. 

A priory, there is no connection between the ET momentum distribution $(p^) of a 
constituent quark model and LF wave function 'il){x^p\). However, the mapping between 
the variables described above turns a normalized solution of the ET equation of motion 
into a normalized solution of the different looking LF equation P3|. Because the ET 
function depends on the relative momentum it is more convenient to use the quark- 
antiquark rest frame instead of the 5-meson rest frame. Recall that in the LF formalism 
these two frames are different. As a result one obtains the LF wave function as ip{C,,p'±) = 
{dpz/d^)^{p']_ +p'j.{C,,p'±))- Exphcit form of this function is given e.g. by Eq. (10) of Ref. 
||T^ . It is wave functions made kinematically relativistic in this fashion, that were used in 



a recent calculation of the Be lifetime [^ . We calculate the photon energy spectra using 

the three representative LF wave functions corresponding to the non-relativistic ISGW2 

25| , ALl pBf, and relativized DSR [^] constituent quark models[|. The main difference 



between the ET wave functions of these models relies in the behaviour at high value of the 



internal momentum, for further discussion see |Q . We believe that the spread of results 
obtained for these distribution functions is a fair representation of model dependence 
resulting from the inclusion of Fermi motion. 

Having specified the non-perturbative aspects of our calculations, we proceed to present 
numerical results for the photon spectrum in the decay B -^ Xg^. In case A we take 



^The ET ISGW2 function corresponds to the Gaussian distribution ^{p^) conventially employed in the 
ACM model with pp = 0.43 GeV. For ALl and DSR models we use simple analytical parametrizations 
of the ET wave functions lESl . 



nib = 4.8 GeV, rrisp = 0.3 GeV and A = 2 as reference values. These values are motivated 
by a study of the b ^ c decays []T^. For the models B to D we use the constituent quark 



masses listed in Table 1. 

The choice of m^ in our approach deserves some comments. The numerical calculations 
using the constituent 6-quark masses show large deviations of the r{B —* X^'y) from the 
free decay rate Fq oc ml {^ 10% for the cases B to D, see Table 2). This signals the 
appearence of the linear l/rrib corrections to the free quark limit. The reason is that the 
constituent quark models usually employ the 6-quark masses that are 300 — 400 MeV 
higher than the pole 6-quark mass. This fact seems to be a subtlety in applying the 
constituent quark model approach to calculate the non-perturbative corrections to the 
B -^ Xs7 inclusive rate. To overcome the uncertainties induced by the constituent mass 
of the 6-quark we use a simple phenomenological receipe that considerably improves the 



situation. Notice that, as in the ACM model |jT9|, l/mj, corrections can be absorbed 

into the definition of the 6-quark mass. We introduce rhb = rub + Srub by imposing the 

1 
condition y{mb) = Xq, where y = J yRLF{y)dy- This condition coincides with that used 



in HQET to define the pole mass of the 6-quark. As a result, the correction to first 
order in l/rrib will be eliminated from the total width in agreement with the general 
statement of HQET. To illustrate our arguments, consider the analytically tractable case 
of the photon spectrum of Eq. (0) with the distribution function given by Eq. (p^SJ). In 



this case /(^) = 2A/(^o-^"i(A)) exp (— A/2[^o/(l — + (1 ~ 0/Co]), and simple analytical 



1 
expressions for Rlf = I RLF{y)dy and y are avaiable, 



Rlf = -(1 - ^0/^2), y = -(1 - 2^0^:2 + ^o^^s), (16) 

Xq Xq 

where k„ = Kn(\)/Ki(\). The 1/Mb correction to Rlf can be absorbed into the 
definition of rhb- Indeed, neglecting ^q and letting rhb ~ (1 — ^0^2)^5, one obtains 
Rlf = 1 + 0(1/M|). The S-meson mass can be eliminated in favour of the b-quark 
mass, so we have the desired result, Rlf = 1 + 0{l/ml). 

We have calculated numerically the values of rhb in different models using Eq. ( p!T]) for 
RhFiy)- Although ifib depend on the assumed shape of distribution, this dependence is 
marginal: the uncertainty on iftb is between 4.6 and 4.7 GeV, depending on the choice 
of f{^iP\)i (see Table 1). These values are consistent whith the h quark pole mass 
^TT'b,poie = 4.8 ± 0.15 GeV ||2^. If we repeat the same exercise by applying Eq. (^ we find 



practically identical values of rhb. 

We first study the photon spectra using Eqs. (H) and ([TI|) . Our results for the photon 
spectra and moments are reported in Fig. 1 and Table 2. The different curves in Fig.l 
correspond to the models A to D. For each case we show separately the spectra calculated 
from Eqs. (|^) and (0) using both rUb and rhb- The influence of the various choices of 
the distribution function and the 6-quark mass can be read off from Table 2, where we 
show both RLF{jnb)i RbF^'nib) and analogous quantities Racm calculated in the ACM 
model, see below. In Table 2 we also show the average photon energy y normalized 
to Mb/'^-i and the moments y'^ and y"^. Here, and henceforth, we present the results 
obtained using Eq. (pi]). The substitution rrib — > rhb modifies the predictions for the 
total rate and the moments y, y"^ by about 8 — 9% for the models B-D, while it affects 
the corresponding quantities by less than 1.5% for the model A. Our final results for the 



integrated photon spectra RLp{fhi,) agree well with the corresponding OPE prediction 

-Rope = 1 + (Ai — 9A2)/2?Tif, where Ai and A2 parameterize the matrix elements of 

the kinetic and chromomagnetic operators, respectively. For Ai = —0.3 ± 0.2 GeV^ and 

A2 = 0.12 ± 0.02 GeV^ -Rope affects the free quark result by a few per cent, -Rqpe = 

0.975 ±0.005. 

The dependence of the energy moments y" on mi, is rather weak, in contrast to that 

1 
of the QCD moments of the shape function < x" >= / x"'F{x)dx, which are very 

-2Mb /A 

sensitive to the difference between Mb and rrib. In particular, changing the 6-quark mass 
from TUb to rhb modifies < x^ > and < x^ > dramatically. Note that the resulting values 
of < x^ > and < x'^ > still have a sizable model dependence. Our predictions for < x^ > 
are somewhat small, although in agreement (for the model D) with the result of Ref. 0] 
and compatible with the results obtained from the QCD sum rules, < x^ >~ 0.5. This 
means that the LF ansatz can be made consistent with the QCD description provided the 
spectator quark is relativistic. This conclusion agrees with that of Ref. |jT9| . 



In order to compare our results with those of the ACM model we have calculated 
the inclusive B —>■ Xgj photon spectra in a simplified ACM model [|19| assuming the 



monochromatic distribution (|T2|) for the free 6-quark. We have used the momentum 
distributions $(p^) of the spectator quark for the models A to D. In all cases we have 
found that the spectra calculated in ACM and LF parton models are almost identical. 
This is not surprising because we have checked numerically that the quark masses m^ 
defined using the L F models practically coincide with the values of the floating b quark 
mass ml = Mb — Jrni + p^ averaged over the distribution <l>(p^). The integrated energy 
spectra Racm for the models A to D are reported in Table 2, they coincide with R^p 
within a per cent accuracy. 

Although we do not consider here perturbative corrections, it is instructive to compare 
theoretical predictions for the Doppler shifted spectrum dB/dEiab in the laboratory frame 
with the CLEO data. To perform a fit to the data, we rebin the boosted photon spectra 
in the same energy intervals as used by CLEO and, for each choice of the distribution 
function, adjust the overall normalization to give the best fit to the data. The results are 
reported in Table 3, and the best fits are displayed in Fig.2. All fits have y^ jridof <^ 1, 
indicating the present accuracy of experiment. Averaging over the models we obtain 

B{B -^ X,7) = (2.5 ± 0.5e..p ± O.S^odei) x 10-^ (17) 

where the last error comes from the model dependence. This result is consistent both with 
the update CLEO measurement |5D| B{B -^ X,7) = (3.15±0.35±0.32±0.26) x 10-^ and 
a recent reanalysis [§ of the CLEO data B{B -^ X,7) = (2.62 ± Ome^ptoil th) x 10"^. 

Finally, we note that the invariant mass Mh of the hadronic final state is related with 
the scaling variable y by Mfj = M^(l— y). Therefore the theoretical results for the photon 
spectrum can be translated into predictions for the hadronic mass spectrum. In Fig. 3 we 
show the invariant mass distribution of the hadrons recoiling against the photon for the 
models A-D. Our predictions for hadronic mass spectra must be understood in the sence 
of quark-hadron duality. The true hadronic mass spectrum for low Mfj ex M|-* may have 
resonance structure that looks rather different from our predictions. A realistic model 
for the hadronic mass spectrum consists of a single peak located at the mass K*{892) 
followed by a continuum (above a threshold value Mt^) which is given by the inclusive 
spectrum and is dual to a large number of overlapping resonances. In Table 3 we show the 



ratios Rk''(892) = B{B ^ K*j)/B{B -^ Xg'j) obtained by the integration of the inclusive 
spectrum in the range Mh < Mth- The result crucially depends on the choice of Mth] we 
use the value of Mth = 1-15 GeV adopted in Ref. |^. Averaging over the different model 
predictions we obtain Rk*{892) = 0.1571q;44. This result agrees both with other theoretical 
predictions and with the CLEO measurement R'^^*(892) ~ 0-17 ± 0.08. 

In conclusion, we have derived a new parton formula, which establishes a simple con- 
nection between the electroweak decay rate r{B —>■ Xg'j) and the rate of a free 6-quark 
decay. The main features of our approach are the treatment of the 6-quark as an on-mass- 
shell particle and the inclusion of the effects arising from the 6-quark transverse motion 
in the B-meson. Our main result is Eqs. (^, (0). Using various 6-quark distribution 
functions we have calculated the photon energy spectra and the corrections to the free 
decay rate. We have shown that the decay width has no linear to l/rub corrections only if 
expressed not in terms of the constituent quark mass but in terms of a mass Thb which is 
defined in the way similar to that used in HQET to define rrib^poie- In this way one avoide 
an otherwise large (and model dependent) correction of order l/rrih but at expence of 
introducing the shift in the constituent quark mass which largely compensates the bound 
state effects. A summary of our results presented in Table 2 shows a fair good agreement 
both with the QCD results and avaiable experimental data. We have also found that the 
photon energy spectra calculated in our LP parton-like approach agree well with the ones 
obtained in the ACM model, provided the same ET distribution function $(p^) is used 
as input in both cases. Finally, we note that it would be interesting to check whether 
the effective values of the 6-quark mass m^ can appear to be approximately the same for 
different channels (6 — > c vs. b ^ u or b ^ s) and for different beauty hadrons. This 
work is in progress, and the results will be reported elsewhere. 
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Table 1 . The values of the constituent quark mass mb and rUsp (in units of GeV) for 
the models A to D. The values of fhb as defined by the HQET condition < x >= are 
also indicated. 



Model 


A 


B 


C 


D 


nib 


4.80 


5.20 


5.227 


5.074 


rfib 


4.73 


4.68 


4.73 


4.60 


Tfisp 


0.30 


0.33 


0.315 


0.221 



Table 2. The total integrated rates Rlf^ Racm, the moments y", n = 1,2,3, and the 
HQET moments < x^ >, < x'^ > of the photon spectra for the cases A to D calculated 



using Eq. (|TT]). 



Model 


A 


B 


C 


D 


RiFirrib) 


0.973 


0.900 


0.905 


0.899 


RLFifhb) 


0.987 


0.986 


0.989 


0.974 


Rkcu{rhb) 


1.008 


1.005 


1.008 


1.003 


y{mb) 


0.873 


0.806 


0.818 


0.795 


vim) 


0.885 


0.873 


0.885 


0.851 


f{mb) 


0.787 


0.726 


0.742 


0.709 


y\ml) 


0.798 


0.778 


0.796 


0.751 


y^irrib) 


0.712 


0.655 


0.676 


0.636 


f{rhb) 


0.722 


0.695 


0.718 


0.669 


< x^ > 


0.378 


0.295 


0.302 


0.456 


< -x^ > 


0.348 


0.130 


0.176 


0.348 
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Table 3. The branching ratios B{B —* ^sl) obtained from the fit to the CLEO data 
and the partial fractions -Ri^*(892)- 



Model 


A 


B 


C 


D 


ro/TB ' 10^ 


2.54 


2.59 


2.50 


2.91 


Rk*{892) 


0.1681 


0.1312 


0.1486 


0.1807 



Figure Captions 

Fig.l. Theoretical predictions for the photon energy spectrum using the LF quark models 
described in the text. Each plot shows the spectra calculated using both the 6-quark 
masses mh (solid lines) and m^ (dashed lines). The thick and thin curves show the results 
obtained using Eq. ([TT]) and (|^), respectively. 

Fig. 2. Theoretical predictions for the photon energy spectrum in the laboratory frame for 
the different distribution functions f{C,,p']_) and the corresponding rhh. The solid, dashed, 
dot-dashed, and dotted curves correspond to the models A to D, respectively. The data 
points show the CLEO results. Both the left-hand plot and the right-hand plot show the 
results of the best fit reported in Table 3. 

Fig. 3. Theoretical predictions for the invariant hadronic mass spectrum for different 
distribution functions /(C,pi). The notations are the same as in Fig. 2. 
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